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£j ; Abstract 

We calculate azimuthal asymmetries and the Callan-Gross .R-ratio for semi-exclusive pion 
production in deep inelastic scattering taking into account higher twist effects. Our results 
■ are qualitatively different from the QCD-improved parton model predictions for semi- 
inclusive deep inelastic scattering. 
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1. The semi-inclusive deep inelastic process 

£+p^e' + h + x, (i) 

where £ and £' are charged leptons and h is a detected hadron, has been recognized long ago 
[1] as an important testing ground for QCD. In particular, measurements of the azimuthal 
angle ip of the observed hadron provide additional information on the underlying QCD 
production mechanism. The angle ip is defined in the target rest frame as follows: The z 
axis is chosen in the direction opposite to the three-momentum transfer q of the leptons, 
and the x-z plane is the lepton scattering plane with the lepton three momenta having 
positive x components. The angle ip is the azimuthal angle of the observed hadron about 
the z axis. 

Different mechanisms to generate azimuthal asymmetries - (cos <p) and (cos 2<p) - have 
been discussed in the literature. Georgi and Politzer [1] found a negative contribution to 
(cos <p) in the QCD-improved parton model and proposed the measurement of this quan- 
tity as a clean test of perturbative QCD. Cahn [2] took into account the effects due to the 
intrinsic transverse momenta of partons bound inside the proton. This nonperturbative 
effect produces a negative contribution to (cosip) and a positive one to (cos2</?). More re- 
cently Chay, Ellis and Stirling [3] combined perturbative and nonperturbative mechanisms 
and were able to fit the measured values for (cos p>) [4] . The nonperturbative effects were 
parameterized in Ref. [3] by Gaussian distributions for the intrinsic momenta of both the 
target (proton) and the observed hadron (pion). In fact, the hadronization effects cannot 
be neglected for the fit of Ref. [3] to the data. The importance of hadronization effects for 
the azimuthal asymmetry was stressed by Konig and Kroll [5] and also by Berger [6] . The 
calculation of Berger is made specifically for the case of single pion production and takes 
into account pion bound-state effects. 

It is an intriguing feature of Berger's higher twist mechanism that it generates az- 
imuthal asymmetries opposite in sign to those of [3,2,1,5]. However, the competing mecha- 
nisms make it difficult to give reliable theoretical predictions for the azimuthal asymmetries 
and their dependence on the whole set of kinematic variables of the process (1). In this 
paper we reconsider Berger's mechanism and discuss ways of isolating its effects from 
those of the competing mechanisms discussed above. Our motivation is twofold: first, this 
approach calculates and incorporates in a nontrivial way higher twist effects by relating 
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them to the distribution amplitude of the produced hadron; second, a similar approach 
for Drell-Yan process was shown in [7] to be of importance in resolving a long standing 
problem in explaining the data for the angular distribution of the produced lepton. 

2. First we discuss the kinematics of process (1). Let / and V be the four-momenta of 
the leptons with q = 1 — V being the four-momentum transfer and Pjv and Ph those of the 
proton and the observed hadron. We define the usual set of kinematic variables 



Q 2 P N -q Ph- P. 



X B = 777^ > y = ~F» r> z = TT~ ' ( 2 ) 



2P N -q' P N -V q-Pi 



N 



where Q 2 = —q 2 . We will also use Q = \[Q 2 and v = q°, E = 1° defined in the 
proton rest frame. The transverse momentum of the observed hadron is given by = 
Pt (cos <p, simp), where ip is the azimuthal angle introduced in the previous section. The 
differential cross section for semi-inclusive deep inelastic scattering is given by 

2Pl da _ 2na 2 My 



d?P h dx B dy Q 4 



L^W", (3) 



where M is the proton mass and and W^ v are the leptonic and hadronic tensors 
respectively. It is convenient to introduce helicity structure functions, 

W T = l(6%W^ x + 6%W^ y ), 
W LT = -(e£W lil ,<Z + e»W llv e v x ), 

WTT = l(e%W^6»-e%W^ y ), (4) 

where e x , e y and are polarization vectors for the vector boson polarized respectively in 
the x, y and z directions. Following Ref. [8] we introduce the following hadronic structure 
functions, 

Hl = M W 

2z 

H * = (W l + Wt), 

2z\q\ 2 



The differential cross section is conveniently expressed in terms of these structure functions 
[8], namely in the deep inelastic limit we get 

Q 2 da Ana 2 ME 



(x B y 2 H x + (1 - y) H 2 



2 

+ ^(2-y)^/T^ H 3 cosy? + ^(l-y)H 4 cos2p), (6) 



dxsdydzdp^dip Q L 

A 

where Hi = Hi(x B , z.p 2 ^, Q 2 ). 

This decomposition exhauts the kinematical information for reaction (1). Apart from 
the overall normalization factor one may specify three independent observables. Two of 
them we take to be the angular asymmetries, 

{OOSLP) = 2 ~Q x B y 2 H[ T {l-y)H,> (7) 

(C ° S2V?) = 2 Q 2 XB y 2 Hl + (l-y)H 2 (8) 
For the third observable we propose 

R = H2 ~ l XB Hl , (9) 
Hi 

which is the analog of the Callan-Gross .R-ratio [9] in inclusive deep inelastic scattering. 

3. We will now calculate observables (7-9) for the specific case of semi- exclusive pion 
production. More precisely, we want to consider a process in which a direct pion with 
z close to 1 observed in charged lepton - proton scattering is not surrounded by other 
high-pj- hadrons. In other words, the pion is not a member of a jet. The events we are 
interested in thus form only a small part of the whole semi-inclusive sample. If none of 
the specifications above apply to the observed hadron it is expected that the usual QCD- 
improved parton model diagrams of Refs. [1,3] would give the dominant contribution. 
Most of the experiments performed [4,10] did not make these specifications*. 



* Note that in the first paper of Ref. [10] single pions were identified and a nonzero asymmetry 
was observed. However, the data were not fitted to the general form of the (/^-dependence given 
in eq. (6). Thus it is difficult to draw conclusions about the asymmetries defined in eqs. (7), (8). 



The dominant contribution to the semi- exclusive pion production described above is 
given by the diagrams of Fig. 1 a,b. Here the amplitude A for the reaction 



u + e — > e + n + d, 



(10) 



is obtained [11,6,7] by convoluting the hard scattering partonic amplitude T(u + e 
e~ + ud + d) with the pion distribution amplitude 0(£), 



This factorized form of the amplitude should be valid for z close to 1. 

For the hadronic cross section (6) we have to convolute the cross section obtained 
from the amplitude of eq. (11) with the parton distribution function of the proton (see 
[6,7] for more details). A typical contribution to the hadronic tensor W^, that is to the 
hadronic structure functions Hi, is depicted in Fig. 1 c. 

We note that no intrinsic transverse momentum for the pion constituents is introduced 
in our model. The source of the observed pr of the pion is the gluon exchange, see Fig 
1. For the sake of simplicity of our presentation the intrinsic transverse momentum of the 
proton constituents will be neglected in the formulae below, but will be discussed later. 

For the structure functions of eq. (5) we find 




(11) 



Hi = M— ({I 2 (z,p T /Q)] 2 + gj[/i(WQ)] 2 ), 
H 2 = M ^ [I 2 (z,p T /Q)] 2 + A^z 2 [h{z,p T /Q)] 2 + ^[h{z,p T /Q)] 2 ), 

2 

H 3 = N2zI 1 {z : p T /Q) (l 2 (z,p T /Q) - pj^pr/Q)), 
H 4 = -N2h{z,p T /Q)h{z,p T /Q), 



(12) 



where I\(z,pt/Q) and I 2 {z,pt /Q) are defined as follows, 



h(z,p T /Q) 



h{z,p T /Q) 




az 2 



z 2 I 1 (z,p T /Q). 



P 2 t/Q 2 ) ' 



(13) 
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Since we will only be interested in the ratios (7)-(9), the overall factor Af is irrelevant. 

Using crossing invariance the results for Hi can also be obtained by analytic continu- 
ation from the formulae for the angular distribution coefficients for the Drell-Yan process 
given in [7]. The transformation formulae are given in the Appendix. 

Substituting the expressions for Hi given in eq. (12) in eqs. (7)-(9) we obtain predic- 
tions for our observables for a specified distribution amplitude </>(£). Since the expression 
for H s turns out to be positive and the one for H4 to be negative, we find (cos (p) > and 
(cos2</?) < 0. 

Now we want to discuss the effects of intrinsic transverse momenta kx of partons 
bound inside the proton, since from the work of Cahn [2] we know that such effects alone 
produce asymmetries opposite in sign to our results. We introduce intrinsic kx by modi- 
fying the distribution function of the proton, 

1 f tv k 2 1 

G q/N (x,Q 2 )^d 2 k T ^j^ exp j-^pjG./^Q 2 ). (14) 

In the diagrams of Fig. 1 the transverse momentum of the quark from the proton is 
now set equal to k T and the quark is off-shell. More precisely we use the parameterization 
Pq=x p~ =0, Pq= k T . The use of different prescriptions will modify our results 
in subleading orders of kr/Q- In the calculation we use the mean value theorem for the 
integration over the absolute value of k T and neglect subleading terms in kr/Q- 

We plot our results as a function of z for three different values of (kx)- In Fig. 2 we 
plot the observables of eqs. (7)-(9) for Q = 2.5 GeV, pt = 0.5 GeV and y — 0.4. For the 
distribution amplitude of the pion we choose a two-humped form, 0(£) = 30(1 — £)£(1 — 
4(1 — £)£) [12]. Here we neglect the evolution of </>(£) with Q. The solid line corresponds 
to (k T ) = 0, the dashed line to (kr) = 0.25 GeV and the dash-dotted line to (kx) = 0.5 
GeV. In Fig. 3 we plot the same observables for Q = 5 GeV and pt = 1 GeV. The choices 
of (kr) and y are the same as in Fig. 2. Figs. 2 and 3 depict our main results. 

An important feature of our model is a steep increase of the i?-ratio as z — > 1. This is 
analogous to the steep decrease of the A parameter in the end point region in the Drell-Yan 
process predicted by the same higher twist model [6,7]. This behavior of A was observed 
experimentally [13]. 

As expected, a nonzero (kr) causes a change in the prediction of our model for the 
azimuthal asymmetries and the Callan- Gross .R-ratio. This change becomes insignificiant 
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for higher values of Q : see Fig. 3. For Q high enough, the analytic expressions (12), (7)-(9) 
can be used even for nonzero (kx). 

In Fig. 4 we show the dependence of the .R-ratio and of the azimuthal asymmetry 
(cost/?) on different distribution amplitudes assuming vanishing (/ct), Pt/Q = 0.3 and 
y = 0.4. To demonstrate that our predictions depend essentially only on the value of 
((1— = Jo 1 an d not on t ne specific shape of the distribution amplitude, 

we choose the two-humped distribution amplitude of Figs. 2, 3 (solid line in Fig. 4) and 
a convex amplitude 0(f) = [(1 - 0£] 1/3 [B(4/3, 4/3)]" 1 (dashed line) which both have 
((1 — f) -1 ) = 5. Here B(x, y) denotes the Beta function. For comparison we plot results 
also for the asymptotic amplitude </>(£) = 6(1 — f)f (dash-dotted line) with ((1 — f) _1 ) = 3. 
The magnitude of the .R-ratio in the endpoint region and the value of (cos ip) are decreasing 
with increasing ((1 — f an< i are no ^ sensitive to the shape of the distribution amplitude. 

As demonstrated, higher twist effects may be isolated in semi-exclusive pion produc- 
tion for moderate values of Q. We therefore look forward to future precision experiments. 
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Appendix 



By using crossing invariance we can relate the hadronic structure functions Hi for semi- 
inclusive deep inelastic scattering defined in eq. (6) with the angular coefficients for the 
Drell-Yan process defined in the Gottfried- Jackson frame by 

2 dcr v 
= Q 9 a - ~ (1 + A cos 2 6> + ,u sin 26> cos ^ + - sin 2 9 cos 20) , (Al) 
dq z dqj,dxdil 2 

where the angular distribution coefficients A, p and v are functions of the kinematic vari- 
ables x = (q ■ Pn)/(P-k ■ Pn) and q\jq 2 . Here, P n and Pn are the momenta of the initial 
state hadrons and q is the four-vector of the produced photon with transverse momentum 
q T - 

We get the following expressions relating coefficients A, p, and v of the Drell-Yan process 
with the structure functions Hi of eq. (5): 

x B H^pt/Q) = -N[l + X(x',p) - A(x',p)}/2 

H 2 (z,p T /Q) = —2 N [\(x',p) - ^A(x',p)} 

VT IQ H 3 (z, Pt /Q) = N [-2^ + AV.p)] 

VP 



with 



and 



p 2 T /Q 2 H A (z, p T /Q) = N [u(x', p) + A(x', p)} ( A2) 



A ( x '' P) = 77~T~^2 [ A ( x '' P) ~ "(x'i P)l 2 ~ ~ A ^ P( x 'i P% 
I 1 + P ) v p 



A'(x', p) = (1+ ^ 2)2 [(1 - p 2 )(\(x', p) - V {x\ p)/2) + 4^p(x', p)] (A3) 

x'(z) = -, p{z,p T /Q) = i -^!Q , (A4) 

z z 



and N is an overall normalization factor. 
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Figure Captions 

Fig. 1: Diagrams (a) and (b) give the leading contribution to the amplitude of reaction (10). 
Diagram (c) gives a typical (one out of four) contribution to the hadronic tensor W^ v . 

Fig. 2: The observables R, (cost/?) and (cos2<p) defined in eqs. (7)-(9) vs. z for Q = 2.5 
GeV, p T = 0.5 GeV, y = 0.4 and 0(f) = 30(1 - f)f (1 - 4(1 - 00- The solid line 
corresponds to (kx) = 0, the dashed line to (kr) = 0.25 GeV and the dash-dotted line 
to (k T ) = 0.5 GeV. 

Fig. 3: The same quantities as in Fig. 2 for Q = 5 GeV and pt = 1 GeV. The choices for (kr) 
and y are the same as in Fig. 2. 

Fig. 4: The .R-ratio and (cost/?) as a function of z for (kr) = 0, pr/Q = 0.3 and y = 0.4. 
The results are plotted for three different distribution amplitudes. The solid line 
corresponds to the two-humped distribution amplitude also used in Figs. 2, 3. The 
dashed line corresponds to the convex amplitude = [(1 - 0^] 1 / 3 [i?(4/3, 4/3)] -1 
and the dash-dotted line to the asymptotic one </>(£) = 6(1 — 
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